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Abstract 

The ac nonlinear dielectric response X3( w > T) of glycerol was measured close to its glass transition 
temperature T g to investigate the prediction that supercooled liquids respond in an increasingly 
non-linear way as the dynamics slows down (as spin-glasses do). We find that x^i^^T) indeed 
displays several non trivial features. It is peaked as a function of the frequency oj and obeys 
scaling as a function of cjt(T), with t(T) the relaxation time of the liquid. The height of the peak, 
proportional to the number of dynamically correlated molecules N corr (T), increases as the system 
becomes glassy, and X3 decays as a power-law of uj over several decades beyond the peak. These 
findings confirm the collective nature of the glassy dynamics and provide the first direct estimate 
of the T dependence of N corr . 
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Most liquids undergo a transition to an amorphous solid state, the glass, when the tem- 
perature T decreases to their glass transition temperature T g . This transition is similar for a 
vast variety of systems such as polymeric, colloidal, molecular liquids. This ubiquity echoes 
the universality of critical phenomena, where the emergence of long range order makes irrel- 
evant most of the microscopic details, yielding the same critical behavior in various systems. 
This is why it has been argued for a long time that some collective effects should be asso- 

n 

dated to the glass transition, ultimately related to the proximity of a phase transition 
An immediate problem with this idea is that standard equilibrium correlation and response 
functions remain those of a normal liquid when T — y T„ [2J|. In the two past decades, how- 

n n 

ever, the heterogeneous nature of the dynamics close to T g was established |3|, |4J, and a 
length scale associated to these dynamic heterogeneities was estimated f|. The increase 
of the number of correlated molecules N corr when T decreases towards T g is expected to 
explain the main aspect of the glass transition, i.e. the huge increase of the relaxation time; 
N CO rr{T) has therefore become a central concern in the field. The first estimate of this de- 
pendence js S] relied on the temperature derivative of a two-point correlation function. Its 
relation with the number of dynamically correlated molecules N corr (T) is however ambigu- 
ous because it involves a temperature dependent prefactor difficult to estimate and control 
. Furthermore, this estimator leads to a low-temperature divergence of N corr (T) for a 
purely Arrhenius system, for which no collective behaviour is expected, at least at first sight 
[71] . An indisputable experimental estimate of the T-dependence of the dynamical correla- 
tion volume is thus still lacking. In this work, we fill this gap by measuring the anomalous 
increase of a physical susceptibility. It leads to a direct estimate of N corr (T), up to a nu- 
merical prefactor that is now temperature independent. Our findings show that supercooled 
liquids respond in an increasingly non-linear way approaching the glass transition. For many 
systems, e.g. spin glasses, such an increase is related to criticality. In the present case, it 
suggests that an underlying phase transition could possibly be present as well, although it 
does not necessarily imply it. 

The spin glass transition taught us that two-point functions can be blind to the amorphous 
long range order and that this critical behavior is revealed in particular by the third order 
nonlinear magnetic susceptibility 10, [ll]. It was recently argued 12, 13] that the frequency 
and temperature dependent nonlinear dielectric susceptibility X3( w ) T) in supercooled liquids 
plays a role similar to the nonlinear magnetic susceptibility in spin glasses and that its 
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increase when T decreases would be a signature of the incipient long range amorphous 
order. Accordingly, X3( w ) T) should display a peak for ojt(T) ~ 1 whose height grows as one 
approaches T g . Here r(T) is the relaxation time of the supercooled liquid which increases 
rapidly as T decreases toward T g . Contrary to spin glasses, however, Xzi 1 ^ — 0, T) should 
remain trivial since at long times glasses still behave as disordered liquids [3] . One expects 
that when r(T) is large, X3 takes the following scaling form: 

XMT) « e °^ 2Q3 iV corr (T)^(u;r), (1) 

where N corr (T) is the T-dependent average number of dynamically correlated molecules, 
Axi = Xi( u = 0) — ~~ °°) is t ne P ar t °f the static linear susceptibility corresponding 
to the slow relaxation process we consider, a 3 the volume occupied by one molecule, and % 
a certain complex scaling function that goes to zero both for small and large arguments. The 
humped shape of \H (cot ) | i s a distinctive feature of the glassy correlations. Indeed, in the 



'no correlation case' 
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15], X3( w > T) is given by the prefactor of Eq. ([TJ times a function 
that, for all T, reaches its maximum value at u — 0. Furthermore N corr (T) does not show 
any temperature dependence in that case. 

We have devised an experiment to measure, for the first time, Xzi^iT) of a classic glass- 
former (glycerol), via third harmonics measurements. Our findings are in good agreement 
with the predictions of Eq. (CTJ: the scaling, the growth of the number of dynamically corre- 
lated molecules with decreasing T, the humped shape of \7i\ and its power- law dependence 
for wr 3> 1 are all confirmed experimentally. The anomalous increase of X3 is compatible 
with the existence of incipient critical fluctuations when T — > T g , even if our measurements 
are too far from the possible critical point to be a clue of criticality. 

Our measurements were performed at temperatures between T 9 + 4K and T 5 + 35 K, with 
T g 190 K. The linear susceptibility = Xii^) + iXii u ) quantifies the first harmonic 

response of polar molecules to a periodic excitation field Ee~ tuJt . The polarization, P(u)e~ luJt , 
of the system is given to first order by P(oo) = e^Xii^E. We define the relaxation time 
t(T) as the inverse of the frequency f a where Xi is maximum. The nonlinear response is 
dominated by Xsi 00 ) which give s to first order the magnitude of the third harmonics response 
to the field at pulsation u [l^. Q, P(3u)e~ 3lu,t = eox^)[Ee' lu)t Y (the second harmonics 
is zero because of the E — > — E symmetry). While xii^iT) has been widely studied in 
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FIG. 1: (Color online) For each of the five temperatures labelling the curves, X(lj, T) = |x3( w > T)\ x 
/csT/[eo(Axi) 2 a 3 ] is given as a function of frequency / = u/2tt. The arrows indicate the five 
relaxation frequencies f a (T), for which Xi(w) is maximum. The dashed lines are guides to the 
eyes. 



supercooled liquids [17H19|. the T and uj dependence of X3 nas n °t been measured so far, 



for any liquid, glassy or not [20 ]. 

In a closed cell, two high purity glycerol samples were prepared between gold plated 
copper electrodes 2 cm in diameter, with mylar® spacers ensuring inter-electrode distances 



of 19 /im and 41 /im 
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The cell was placed in a cryogenerator, and the temperature T 
of the samples was regulated with a precision of 50 mK. A low harmonic distortion voltage 
source yielded a field E ~ 7 x 10 5 V/m (rms) in the thinnest capacitor, and x 3 (o;,T) was 
obtained from the measured third harmonics current I(3u). This nonlinear signal was 10~ 7 
to 10~ 5 of the linear one, well below the harmonic distortion of electronic devices. To get 
rid of them, we used a high sensitivity method based on a bridge containing the two glycerol 
capacitors [16] . 
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FIG. 2: (Color online) The quantity |%3(a;,T)| at five temperatures (same data and symbols as in 
Fig. 1), normalized to its maximum value at each T, plotted as a function of the ratio / /f a (T). 
The straight line is a power-law with an exponent -0.65. Inset: corresponding phase of X3(w,T). 



Figure 1 shows, for five temperatures, the u dependence of X(u,T) = ^(u;, T)| x 
/eg? 1 / [(Axi) 2 a 3 e ]. According to Eq. (pQ), this quantity should be N corr {T) \H(lot)\. Two 
striking results, confirming the predictions of Eq. (JT]), can be seen in Fig. 1. First, X(u,T) 
is actually peaked at a frequency /* of the order of l/r(T); more precisely /* ~ 0.21/r(T). 
Second, the height of the peak, which should be proportional to N corr (T), increases signif- 
icantly as T decreases. The small value of X(u,T) can be understood from calculations of 

nn 

X3(a>,T) for molecules undergoing independent rotational Brownian motion [14|, [15[, which 
should be accurate for simple supercooled liquids at high T. They indicate that for the 'no 
correlation' case, X(cj, T) < 0.2. 

To what extent is the scaling predicted by Eq. (TjQ) verified by x 3 (u;,T)? Fig. 2 shows 
the same data as in Fig. 1 but plotted as a function of log(/// a ) with f a = r(T) -1 , and 
normalized by their maximum at each T. We also show the phase of xs- Both the modulii 
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FIG. 3: (Color online) The measured quantity max w [X(u), T)] = max w (|x3(w, T)\) x 
A;_BT/[eo(Axi) 2 a 3 ] plotted vs. temperature. The continuous line is the number of correlated 
molecules estimated from T\T (see Q) with an arbitrary normalization chosen so that it coin- 
cides with the experimental points at 202 K. 

and the phases collapse fairly well on a single master curve, as predicted by Eq. ([T]). A weak 
departure from scaling occurs at low u for the highest T. This is not surprising: scaling 
should not hold far above T g , where the dynamical correlations become short-ranged 12L Il3|. 
urthermore, the nonlinear response results both from trivial dielectric saturation effects 
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22] (always present for lot < 1), and from the non trivial dynamical correlations 
contribution. As the latter vanishes when oj — >■ (jl^, [3]), the contribution of the trivial 
saturation effect should dominate at high T and low oj. The observed departure from scaling 
can thus be explained by the different T and u dependencies of the two contributions. 

Figure 3 gives the T dependence of the maximum value of X(u, T), reached for u = u>* = 
27rf*. Since scaling is well obeyed, Eq. ([1]) tells us that this maximum value is proportional 
to N corr (T). A clear increase of this quantity as T decreases toward T g is visible in Fig. 



3. This is one of the main results of our study, and a strong experimental evidence of 
a growth of the dynamical correlation length close to T g . We now investigate how our 
results are related to another estimate of the size of the dynamical heterogeneities, using 
the temperature derivative of Xi( u )- 

I n |z, S], the T-derivative of the two-body correlation C(t) (or of was introduced 

as a new response function xr(t) = dC{t)/dT. It turns out that Txrif) is also propor- 
tional to a dynamical correlation volume, but with an unknown, and possibly T dependent 
prefactor [3, Q]- The line in Fig. 3 is the number of correlated molecules estimated from 
max^(T(9(x / 1 (a;)/Ax 1 )/9T) where, due to the unknown value of the normalization was 
chosen such that it coincides with X(u*,T = 202K). We see that the T dependencies of 
these two quantities are close. However, whereas the increase of X(u*,T) is a proof of a 
growing correlation length, one needs some extra assumptions to infer this from xt Q] ■ That 
X3 and xt lead to a similar T dependence for N corr (T) suggests that these assumptions are 
indeed warranted (at least for glycerol) and validates the simpler procedure advocated in j(| , 
and used extensively in [8(, to extract N corr (T). 

A precise identity relating xt and X3 nas been obtained in [l^]. It holds at low enough 
oo when the linear response depends only on external parameters (T, density, electric field, 
. . . ) via the dependence of t(T) on these parameters, a property called time temperature 
superposition (TTS), and reads: 

X 3 (w) ~ kx t (2u) ut < 1 , (2) 

where k is independent of u and xt is the Fourier transform of xt- Corrections to this 
relation in the regime wr <C 1 are due to weak violations of TTS and to the additional term 
that corresponds to u independent non-cooperative degrees of freedom. The continuous line 
on Fig. 4 presents X(u,T) calculated from the linear susceptibility, according to Eq. ([2]) 
with k = 7.6 x 10~ 17 K(m/V) 2 . This value of k was chosen to be consistent with Fig. 3, 
i.e. by demanding that the maximum over u of k |x^(2o;)| corresponds to X = 0.40 which 
is the value of the continuous line of Fig. 3 for T = 210. 2K. Although Eq. ([2]) should 
only be valid at small u, we see that it reproduces qualitatively the overall behavior of the 
data at 210.2 K. At low frequency, ojt < 0.2, the discrepancy can be understood by the 
contribution of the u independent trivial saturation effect noted above. For lot > 1, our 
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FIG. 4: (Color online) The measured quantity X(u, T) at 210.2 K (filled squares with dashed line) is 
compared to the same quantity calculated by using Eq. (0) with k ~ 7.6 x 10~ 17 K(m/V) 2 ). Dotted 
line: Calculated heating contribution (upper limit, the arrows indicate the possible reduction [bH]). 
Inset: The phase of X3(w, T) at 210.2 K (same symbols as for the main figure, dotted line not 
shown) . 

data exhibit a clear power-law behavior with a fitted exponent —0.65 ± 0.04, see also Fig. 
2, which is close to that describing the decay of Xr(2w). This behavior is predicted by 
Mode-Coupling Theory [l3], and our experimental exponent is compatible with the MCT 
prediction —b ~ —0.6 for glycerol j^j. Note however that MCT should in fact be relevant 
only at temperatures much higher than the ones we focused on. Hence, the existence of a 
power-law regime for large u with identical exponents for both xi an d X3 appears to be a 
generic property, valid outside the MCT regime. We think that this might be related to the 
possible incipient criticality of the system. 

We now address the question of a possible heating contribution to our results. Such ef- 
fects would come from the fact that applying the field E(t) = Ecos(ut) across the glycerol 
sample leads to a dissipated power per unit volume p(t) oc uXi(oo)E 2 (l + cos(2u;i — (f>)) 
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sec 



14j). The resulting temperature increase is the sum of a constant term and a term 



8T2u(t), oscillating at 2ou. This 5T2 W (t) oscillation leads to a modulation of the linear sus- 
ceptibility 5xi(u)(t) = (dxi(u))/dT)5T 2uj (t), thus to a 3u modulation of the polarization 
P(t) oc XiE- We calculated precisely this 3a; heating contribution using the heat propaga- 
tion equation [ijj]. This calculation gives an upper limit to the heating effects contribution, 
in particular for / > f a , because the finite relaxation time r of the molecules prevents them 
from following instantaneously the temperature oscillation 5T2 LJ (t) (as already advocated by 



Richert et al. 



24|). A representative example of these calculations is given on Fig. 4 : the 



heating effects are negligible in most of the frequency range of interest, in particular around 
the Xzi^) peak. Thus, our estimate of N corr (T) and the power-law dependence of Xzi 1 -^) are 
not affected by heating. 

Our results should allow significant progress of models which predict or assume dynamical 
correlations in supercooled liquids. Such models should not only reproduce the increase of 
the correlated volume we found as T — > T g (see Fig. 3), but they should also account for 
the magnitude and shape of X3(^,T), given by the H, function in Eq. (TB, which carries 



an important information on the physics of the dynamica 



corre 



ations 



13| . Recently, 

Richert et al. put forward a phenomenological model [21|, [24J, [25( which accounts for their 
nonlinear susceptibility data at lu. In this model, the supercooled liquid is thought of as 
a collection of "Debye-like" dynamical heterogeneities (dh), each of them having its own 
relaxation time r^. The model posits that the electrical power absorbed by each DH raises 
its temperature above that of the phonon bath. A non-trivial and crucial assumption is 
that heat exchange between the slow degrees of freedom and the phonon bath is set by the 
local relaxation time t^h, and not by a microscopic vibration time. As a result, each DH 
has its own Active temperature with a dc and an ac component. Again, the ac component 
5T dh ^uj{t) leads to a 3u "heating" contribution to the polarization P(t). Such a contribution 
to X3 should however not be considered to be in competition with the critical dynamical 
correlations related to the glassy dynamics [3, [3]. As this model assumes a priori the 
existence of dynamical heterogeneities, it should be seen as a phenomenological description 
of the influence of dynamical correlations on the nonlinear susceptibility. We calculated 
the prediction of this model for the 3a; nonlinear response and found that ^3(0;) is indeed 
peaked at a frequency of the order of l/r(T). However, some assumptions are needed to 
generalize the model to the 3a; response. These assumptions influence the magnitude and 
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position of the peak and deserve further scrutiny. A detailed comparison of this model 
to our data will be presented in a subsequent paper. In any case, we believe that for a 
proper prediction of Xzd^, T) close to T g , a theory of supercooled liquids able to account for 
dynamical correlations is required. 

To conclude, we have provided the first direct experimental evidence that a supercooled 
liquid responds in an increasingly non-linear way approaching the glass transition. By 
measuring the frequency dependent third harmonics response X3(cu,T) to a periodic electric 
field, which is tightly related to the dynamical correlation length, we showed that the number 
of correlated molecules increases as T decreases towards T g , confirmed the validity of previous 
estimates, and found that X3 scales as a function of lot. This opens a new path for probing the 
spatial correlations in both fragile and strong supercooled liquids and in the aging regime of 
glasses and spin-glasses, by systematic studies of non-linear responses. Future investigations 
along these lines might help to unveil the possible critical nature of the glass transition. 
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Evidence of growing spatial correlations at the glass transition 
from nonlinear response experiments 



For the sake of completness, we detail in the four sections hereafter: (i) the connexion 
between the third harmonics of the fundamental response and the nonlinear part of the 
response; (ii) what can be expected for the third harmonics response of a dielectric liquid 
where no (glassy) correlations are present; (iii) some additionnal information about our 
experimental setup; (iv) how the upper limit of the heating contribution was computed 
(Fig. 4 of the letter) and the reduction which can be expected from a more refined estimate. 

I) NONLINEAR RESPONSE MEASUREMENTS THROUGH THIRD HARMON- 
ICS DETECTION 

As in the main paper, we consider the nonlinear response of a dielectric system to a time 
dependent electric field E(t). The most general relationship between the polarisation P(t) 
and the excitation E(t) is a series expansion in E pj. For a purely ac field, the even terms 
are forbidden because of the E(t) — > —E(t) symmetry, which yields : 



In this equation eo is the dielectric constant of vacuum, xi the linear susceptibility and % 3 
the cubic nonlinear susceptibility in the time domain. The dots in Eq. [3] indicate an infinite 
sum involving higher order nonlinear susceptibilities x^i etc- Note that causality implies 
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X3 (t-t' 1) t-4,t-t' 3 )x J E;(t / 1 )£;(4)£;(4)dt , 1 «+.... (3) 
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Xi(t < 0) = 0. The Fourier transform of Eq. [3] for a purely ac field E(t) = Ecos(ut) gives: 



eo 



+ 



E 
~2 
E 
2~ 

E 



3E 2 

Xi{u) + -r-X3(-w,a;,a;) + ... 
3E 2 



5(u/ — w) 



^(a;' + u;) 



E 3 



(4) 



where the polarization P and the susceptibilities x« are now taken in the frequency domain 
and the dots indicate again infinite sums involving higher order terms. The response P(t) 
to E(t) = E cos(ut) can thus be written 

P(t)/e = Re [(Exi{«>) + 3/4£ 3 X3 M + ...)e~ 1ut ] + Re [l/ AE^x^e^ 1 + ...] + .... (5) 

To obtain Eq. (jSJ), we have used the fact that because xi an d X3 are rea l i n the time domain, 
their Fourier transform verify Xi(^) = Xi(~ u ) an d xtiwi, ^2, ^3) = X3(~ w i) —^2, —^3) (the 
star denotes the complex conjugate), and the invariance of X3 by permutation of its argu- 
ments. For simplicity, we write X3( u ) — Xsi^^^) an d = X3(~ uj,u,u). Eq. (j5J) can 
be written 



P(t)/e = E(x'i cosut + Xi sin cut) + 3/4E (x'3 cos cut + x'i sin cut) + ... 

+ l/4£ 3 (x3Cos3wt + X3sin3u<) + (6) 

where the susceptibilities Xi are given as a function of their real and imaginary parts x\ 
and x'l- F° r practical applications, the modulii and arguments \x%\ a nd Si are rather used: 



P(t)/e = E cos(wt - Si) + 3/AE 3 \ X s\ cos(ut - S~ 3 ) + 



+ ... + 1/AE 3 |xs| cos(3wt - S 3 ) + ... (7) 



We see in the rhs in Eqs [5H3 that the nonlinear susceptibility of interest, namely X3 



13 



appearing in Eq. 1 of the letter, is directly given by the measurement of the third harmonics 
of the polarisation. 



II) THIRD HARMONICS FOR A LIQUID WITHOUT GLASSY CORRELATIONS. 

We briefly summarize here what can be expected for the third harmonics in the case o " 
a liquid without any glassy correlations. We use the work of Dejardin and Kalmykov {2] 
which studies the nonlinear dielectric relaxation of an assembly of rigid polar symmetric 
particles. The particles are assumed to be diluted enough in a non polar solvent, so that 
any interaction between them can be safely neglected. A strong electric field is applied to 
the system, which biases the (noninertial) rotational Brownian motion of the polar particles 
and yields a net dielectric response. This response is expressed by expanding the relaxation 
functions as a Fourier series in the time domain. Then the infinite hierarchy of recurrence 
equations for the Fourier components is obtained in terms of a matrix continued fraction. 

We focus on the case where the field is purely ac, namely E(t) = E cos ut. In the following 
equations, the subscript "D" refers to the calculation of Dejardin and Kalmykov. By using 
their Eq. (21), one gets : 

X3,£>M = 3 - 17u 2 T% + iOJT D {U - 6U 2 T 2 D ) 

X3 , D ^ = 0) (1 + co 2 r 2 D )(9 + 4u,V|,)(l + 9u,V£) ' 1 } 

where td is the usual Debye time since one finds Xid 00 if^p > i- e - Xi 1S maximum for the 
frequency fn = l/(2nr D ). By using Eqs. (7) and (19) of Dejardin and Kalmykov, one can 
calculate all the prefactors and obtain X3,d(w = 0): 

( a\ e (Axi) 2 l , v 

X3 > d{u = 0) = ^TiV^S ' (9) 

where N is the volumic concentration of polar particles, and where, as in the letter, A^i = 
x'id(0) — Xid(°°)- To get a benchmark for the behavior of Xz,triviai for the 'no glassy 
correlations' case, we set fo = fa as well as A^o = 1/a 3 with a 3 the molecular volume in 
glycerol. We thus obtain from the above equations: 

_ e (Axi) 2 a 3 fl X 3 ,d(w) \ 

X3,trivial{U — — — I -; — - . [W 

k B T \5x3,d{u = 0)J 
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Fig. S 1: The quantity \X trivia i \ = \ Xi ' tr f! al{ \} k ^ T \ is plotted as a function of the ratio ///d- Inset: 
same graph in log-log plot showing that X trivia i decreases as /~ 3 when / ^> fr>. 

The last term in the rhs of the latter equation depends only on / / f D , i.e. f / f a . Thus, 
the latter equation and Eq. (1) of the letter are formally similar: The quantity X trivia i = 
X3 eo™&xrfa^ T P^ a y s ' ^ or ^ e ' no glassy correlations' case, the same role as X = ^T^^fs in 
the main letter. The behavior of X tr iviai is shown on Figs S1-S2. Three main features are 
noticeable: 

(i) On Fig. SI, one sees that \X tr i V iai\ has no peak in frequency, at variance with what 
is observed on Figs (1),(2),(4) of the main letter. This strongly support that, as claimed in 
the letter, the peak in frequency of \X\ is directly linked to the effect of glassy correlations. 

(ii) From Eqs. (!H!)- (|T0l) . one sees that the magnitude of \X trivial \ does not depend on 
the temperature T. The T dependence of X trivia i is only through td(T), since we have set 
fD — fa- This absence of temperature dependence of \X tr i V i a i\ strongly supports our claim 
that the increase of the maximum value of \X\ seen on Figs (1),(3) of the letter comes from 
the variation of the correlation volume when the temperature is decreased towards T g . 

(iii) On Fig. SI, one sees that \X trivia i\ is maximum for u = 0. The trivial part of X3 thus 
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Fig. S 2: The phase of X trivia i is plotted as a function of the ratio ///d- 

dominates over the singular part in the limit of low frequencies / -C f a . Of course the role of 
the trivial part of X3 is larger when the non trivial part is lower, i.e. when the temperature 
is high. This naturally accounts for the departure to the scaling of X as a function of f/f a 
when both the temperature is high and the frequency / is much lower than f a . 

Ill) MORE INFORMATION ABOUT OUR EXPERIMENTAL SETUP 

The principle of the experiment was that of a bridge including two capacitors in which 
the dielectric layer was the supercooled liquid under study (here, glycerol) as described in 
section V of ref. These two glycerol samples with different inter-electrode distances were 
placed in a closed cell filled with high purity glycerol (purity > 99.6 %) in a dessicated 
atmosphere to prevent water absorption. The temperature of the cell was regulated by a 
LakeShore® temperature controller with a 50 Ohms heater fixed on the cell. The cell and 
its heater were placed in a cryostat in which the cold stage remained at ~ 20 K and was 
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connected to the cell through a thermal impedance. This device allowed a temperature 
regulation with a precision of ~ 50 mK, and a temperature quench from 300 to 200 K in ~ 
1.5 hour. 

The electrodes of the two capacitors are 6 mm thick rods of high purity copper on which 
200 nm of gold was evaporated to prevent the formation of any oxydized layer. A highly 
specialised surface treatment reduced any surface defects: The final roughness was about 10 
nm, while the planeity, i.e. the difference between the highest point of the electrodes and 
the mean plane, was better than 1 fim. 

Most of the results presented in the letter were obtained with Mylar® spacers depicted 
in the letter, giving inter-electrode distances of 19 /xm and 41 /^m for the two samples. For 
completeness, we have added a few results obtained very recently by using spacers made of 
resin patterned by photolithography (with a thickness L — 12 /iin for the thin sample and 
25 /im for the thick sample). 

The bridge was fed by a low distortion voltage source (Stanford Research Systems® 
DS360) which gave a rms excitation voltage < 14V at a frequency / = u/2tt ranging from 
10 niHz to 200 kHz. The nonlinear susceptibility X3 was determined after having carefully 
checked the cubic dependence of the measured signal at 3a; with respect to the excitation 
at lu. Last, for a significant fraction of the frequencies and temperatures investigated here, 
we have checked the quantitative consistency between the 3a; signal measured in our "two 



sample bridge" and the one found with the "twin T notch filter" depicted in Ref. 
allows to measure xs by a different method. 



l|, which 



IV) HEATING EFFECTS CALCULATIONS 

Let us consider a sample made of a supercooled liquid excited by an oscillating field 
Ecos(u)t). For small enough E values, the resulting polarisation Pu n is linear and reads : 

^ = X / 1 cosM + X ' 1 'sin(a;t). (11) 

When E is increased, a part of the nonlinear response comes from the dissipated electrical 
power, the volumic density of which is p(t) =Jk<e§)<![ujE' 1 (1 + cos(2o;t — 0)), with = — n + 
2 arctan(x'i/Vi) (see e.g. Appendix of Ref. js]). The resulting heat propagates towards a 



17 



"thermostat" which in our case can be considered to be the copper electrodes (thickness 6 
mm, diameter 20 mm). The resulting average sample temperature increase 5T(t) is 



5T(t) = 5T + 5T 2 cos(2wi - 2 ), 



(12) 



where the mean dc temperature increase STq is larger or equal than the ac one 5T 2 , thus at 
any time ST(t) > 0. <p 2 is a phase shift related to heat transport. As our measurements give 
the nonlinear dielectric response averaged over the sample volume, ST(t) in equation ( jT2j) 
is the temperature increase averaged over the same volume. Using equations ( IH|) - (fT21) . we 
thus obtain for the nonlinear part of the susceptibility due to heating effects: 



P(t) ~ Pun{t) 

e E 



9x,1 6T(t)] cos(a;t) + (^8T(t) ) sin(ut). 



dT 



dT 



(13) 



This is an upper limit on the contribution of heating to the nonlinear response, since we have 
assumed that 5T(t) induces instantaneously a modification of the susceptibility. As already 
advocated in Ref. jsj, this is questionable, specially in what concerns the contribution of 
5T 2 (t) which should be averaged out because of the finite relaxation time r of the system 
(see after eq. [17]). By using equations (I12H13I) . one gets : 



P(t) ~ Plinjt) 

e E 



+ 
+ 
+ 



5T + -5T 2 cos((p2, 



ST, 



o 



-5T 2 cos(4> 2 ] 



dx\ 1 dx'{ 

— + -(5T 2 sin(0 2 )^T cos(wt) 

dy'i 1 <9x'i 
— + ~8T 2 sm(4>2) — sin(wt) 



1 srrp (. ^ d Xl 

-<5T 2 cos(0 2 ) — 



1 dx'l 
-5T 2 sin(</) 2 )— — 

2 y 1 dT 



T d\\ 1 dx'{ 

-5T 2 sin(0 2 )— + -8T 2 cos(4> 2 ) — 



cos(3ojt) 
sin(3u;t). 



(14) 



The four terms in the right hand side of equation (1141) give the nonlinear response of the 
system due to heating that we define in general by 

— eo ff /m< ^ = ( 6 x'i)hCos(ut) + (5xi)hsin(ut) + ^Xz,h cos(3wt) + ^~x'{h sin(3wt). (15) 



In the right hand side of equation ( 1T51) each of the four terms are proportional to E 2 . x'zh 
and x'ih d° n °t depend on E, while (8x'i)h and (SxT)h are proportional to E 2 . Our notations 
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Mi. 



have been chosen in order to remain as close as possible to the notations of refs 

We now calculate the expression of 8T(t) that one has to introduce in equation (fT4l) in 
order to obtain the nonlinear response in equation (flo]) . We just consider the simple heating 
effect, when the dynamical heterogeneities are not considered: The supercooled liquid is 
thus just characterized by its thermal conductivity Kfh an d its total specific heat c. As in 
ref. [6], we consider that c is frequency dependent due to the fact that the slow degrees of 
freedom cannot contribute to c for frequencies much larger than f a = 1/r. Let us define 
(x, y) as the plane of our copper electrodes, with z = for the lower electrode and z = L for 
the upper one. Due to their very high thermal conductivity and to their large thickness (6 
mm), the two electrodes can be considered, to a very good approximation, as a thermostat. 
The temperature increase S9(x, y, z, t) of the supercooled liquid at point (x, y, z) and time t 
thus vanishes for z = and z = L. As the diameter D = 2 cm of the electrodes is typically 
one thousand times larger than L, we may consider that 68 does not depend on (x,y). We 
obtain S9(z,t) by solving the heat propagation equation: 

c— — = Kth —^- + P (t). (16) 

By averaging spatially the solution of equation (JIB"]) , we obtain the 5T(t) to be used in 
equation (1T4T) : 

Oln = 

24/ch, 

at u\ at cos(2u;t - 2 ) 

sm = * V + (W (17) 

where r t h = cL 2 / (K t h^ 2 ) and 02 = <fi + arctan(2o;rth). By using these expressions in equa- 
tion (fbil) and identifying it with equation ffl5l) . we obtain the upper limit to the heating 
contribution to X3 (dotted line in the Fig. 4 of the Letter). 

We now move to the problem evoked above, namely the fact that the finite relaxation 
time t of the dipoles which contribute to the dielectric susceptibility should average out 
the modification of this susceptibility due to the oscillation <5T 2 (t), specially in the case 
wr > 1. As a consequence, the heating contribution to %3 should be reduced by a factor 
R{ut) < 1. For a precise calculation of R(ujt), one should replace Eq. (fl3l) by a microscopic 
equation accounting for the dynamics of the dipoles in the case of a thermal bath where the 
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temperature has an oscillating component, which is of great complexity. For a first estimate, 
we make two very simplifying assumptions: 

(i) We assume that the dipoles have a Debye dynamics with a given characteristic time 
t(T). This is a simplifying assumption in the sense that when T is close to T g , it is well known 
that is "stretched" with respect to a simple Debye law. In fact the Debye dynamics 
holds only at much higher temperatures, where the molecular motions are independent of 
each other, which allows to describe the non inertial rotational Brownian motion by the 
Smoluchowski equation for theprobability distribution function of the orientations of the 
dipoles in configuration space (2J, |7| . After an ensemble averaging of this equation, one gets 
the well known Debye equation for the dynamics of the average polarisation j^j : 

OP 

r— + P = e A X iE cos(a;t) (18) 

(ii) we assume that the main effect of ST 2 (t) is to modulate in time the value of r 
while leaving unchanged the (Debye) dynamics. This can be justified by the fact that 
the temperature oscillation modulates the viscosity rj, thus the relaxation time r which is 
proportional to rj. Considering the heatings of eq. (fT7|) . r(t) is now given by 

r(t) = run + (^) m + (^p) 5T 2 (t) , (19) 

where Tu n is the value of r at zero field. In the following, 5t 2 will denote the amplitude 
of the 2u modulation of r due to ST 2 (t) and corresponding to the last term of Eq. ffl9l) . 
Of course, using Eq. ( 1T9"|) for r(t) assumes that 5T 2 (t) instantaneously fully affects r. A 
thorough modelization of this problem could lead to a more involved expression where 5t 2 
should be weaker than in the above expression. As we shall find that the third harmonics 
is proportionnal to 5T 2 /Tu n , see below Eq. (l2~Tj) . we are lead to the conclusion that our new 
estimate should, again, be slightly overstimated. 

We now insert r(t) in Eq. (TlBl and set : 

P(t) = P Un cos(wt - * Kn ) + 8P X cos(ut - ^i) + P 3 cos(3wt - $ 3 ) + ... (20) 

where Pi in , ^u n , SPi, ^1, P3, ^3 are to be determined. As we are only interested in the onset 
of nonlinear effects, Pu n oc E is much larger than 5 Pi oc E 3 and than P 3 oc E 3 . This allows 
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to neglect higher order harmonics (denoted by the dots in Eq. (120]) ) and to resolve Eq. ( !T8|) 
by identification of the terms which have the same frequency and the same power of E. This 
yields: 

where 6t 2 denotes the amplitude of the 2u modulation of r arising from 5T 2 (t) in the last 
term of Eq. ffl9l) : thus 5t 2 oc E 2 , see Eq. f fTT]) . which yields the expected P3 oc E 3 . 

We now have to compare this result to that obtained if we start from Eq. ( fl3|) and 
use a Debye linear susceptibility. A straightforward calculation shows that in that case the 
modulus of the third harmonics of the polarisation is given by the expression of Eq. ( 12 ip 
divided by a function R(ut) given by : 



%) = ^ + B (22) 
A/1 + (3wr K . ^ 2 



tin J 



As expected -R(w = 0) = 1 and R(ut ^> 1) is smaller than 1. We note that R(ut — > 00) = 
1/3, which comes from the fact that r(t) enters in Eq. (fl8|) as a factor of dP/dt: this gives a 
weight 3ut contrarily to the case where one starts from Eq. f fl3|) where this weight is simply 
ujt. This reduction of the effect of ST 2 (t) on the polarisation can be seen as a first estimate 
of the fact that the dipoles average out the temperature oscillations. Of course, one could 
build a much more thorough model of this effect, but the reduction given by Eq. ( 1221) is 
enough to ensure that the power law regime reported in the main letter is not affected by 
the heating contribution calculated here. Finaly, the lengths of the arrows of the Fig. 4 of 
the main letter are those predicted by Eq. (122]) . 
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